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ABSTRACT 

We have studied the X-ray properties of ageing historical core-collapse supernovae in 
nearby galaxies, using archival data from Chandra, XMM-Newton and Swift. We found 
possible evidence of a young X-ray pulsar in SN 1968D and in few other sources, but 
none more luminous than ^ a few 10'^^ erg s~^. We compared the observational limits 
to the X-ray pulsar luminosity distribution with the results of Monte Carlo simulations 
for a range of birth parameters. We conclude that a pulsar population dominated by 
periods ^ 40 ms at birth is ruled out by the data. 
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1 X-RAY SUPERNOVAE AND PULSARS 

X-ray imaging missions have detected ~ 40 young core- 

■ collapse supernovae (SNe); this has led to a deeper under- 

■ standing of the initial phases of shock expansion. However, 
only ~ 10 SNe have been detected in the X-ray band at ages 
> 20 yr, even though there is a reliable optical record of his- 
torical SNe in nearby galaxies over « 100 yrs. On the other 

, hand, core-collapse SN remnants are found with much older 
' ages 350 yrs). But we still know little about the shock 
[ propagation into the ambient medium during the transition 
. between the SN and SN remnant phases (ages of ~ 30-300 
' yrs). Thus, it is important to search for and study historical 
, SNe as far as possible into this age range. 
' From the X-ray point of view, the main phases of SN 
' evolution are (Fransson et al. 1996): 

a) SN phase — circumstellar medium interaction: the shock 
propagates into the stellar wind of the progenitor star. Emis- 
sion from the forward shock (gas temperature a few keV) 
dominates only in the very early epoch [t ^ 100 d); then, 
after the shell of ejecta becomes optically thin, the reverse 
shock provides the X-ray photons (thermal spectrum at 
kT fti 0.5-1 keV). The older the SN, the further back in 
time we probe the mass-loss rate of the progenitor star, and 
hence its evolution in the last few thousand years. 

b) SN Remnant phase = interstellar medium interaction: 
the expanding shock reaches the boundary between the stel- 
lar wind bubble and the surrounding interstellar medium. 
When that happens, the luminosity decay flattens out, and 
the remnant may keep a constant or very slowly declining 
luminosity and temperature for ~ a few 10"^ yrs. 

This scenario can have an additional component. A 
large fraction of core-collapse SNe will leave behind a spin- 
ning, isolated neutron star (NS) powered by electromagnetic 



losses. For a dipole magnetic field B, the energy loss is given 
by (assuming an orthogonal rotator) Erot = fl^ /6c^ , 
where R is the NS radius and Q its angular velocity. There 
is an empirical relation (Possenti 2002) between the X- 
ray luminosity of the pulsar and its rotational energy loss: 
log Lx, [2-10] = 1.34 log -Erot — 15.34. Analogous relations 
were inferred by Li et al. (2008) and Kargaltsev & Pavlov 
(2008). The X-ray luminosity includes both the magneto- 
spheric emission from the pulsar itself, and the synchrotron 
emission from the pulsar wind nebula. It has a power-law 
spectrum with photon index T « 2. It declines as Lx = 
Lx,o(l + t/to)~'^ , where to ~ a few 10'^ yrs. For t i^, to the 
pulsar flux does not vary signiflcantly. Therefore, for some 
SNe with very energetic pulsars (e.g., the Crab), the hard 
X-ray emission from the pulsar itself becomes progressively 
more important than the thermal emission from shocked gas, 
as the SN ages. In external galaxies. X-ray imaging cannot 
resolve the shocked gas from the point-like pulsar, but X-ray 
spectra and colors can tell which component dominates. 

The pulsar luminosity distribution is determined by 
the magnetic fleld and initial spin period distributions. 
The magnetic fleld distribution is relatively well constrained 
(~ lO'^'^-lO'^'^ Gauss). Instead, the birth spin is very uncer- 
tain. Some radio studies have suggested a large population 
of millisecond pulsars (Arzoumanian et al. 2002), while oth- 
ers have claimed millisecond pulsars are very rare (Faucher- 
Giguere & Kaspi 2006). We have recently proposed a new, 
independent method to constrain the initial period distribu- 
tion of pulsars and solve this controversy (Perna et al. 2008). 
Using Chandra, Swift and XMM-Newton data, we have com- 
pared the observed X-ray luminosities (or upper limits) of 
young SNe with the predicted distribution of pulsar X-ray 
luminosities. 
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Figure 1. The dashed line shows the distribution of 2—10 keV lu- 
minosities (measurements or upper limits) for historical SNe with 
ages > 30 yrs; they are treated as upper limits for the underlying 
pulsar luminosities. The solid line (with shaded Icr error) shows 
the predicted pulsar luminosity distribution. 

2 THE OLDEST X-RAY SUPERNOVAE 

A spin-ofT of our observational study (Soria & Perna 2008) 
was the discovery or recovery of three pre-1970 X-ray SNe, 
that is in the characteristic age range when the outgoing 
shock reaches the interstellar medium. The oldest X-ray SN 
we have discovered is SN 1941C, with Lx ~ 5 x 10^'^ erg s~^ 
in the 0.3-10 keV band (Table 1). For SN 1959D, we infer 
Lx ~ a few lO^'' erg s"\ For SN 1968D, the emitted lumi- 
nosity Lx ~ 2 X 10"^* erg s~^, mostly thermal (i.e., from the 
shocked gas), but with an additional power-law component 
with photon index « 2 and Lx « 2 x 10^^ erg s~^. In total, 
there are only 27 historical core-collapse SNe from 1900- 
1970, at distances ^ 15 Mpc (for many other old SNe, the 
type is undetermined). Of them, 17 are so far undetected 
in X-rays, down to ~ a few 10"^^ erg s~^; 4 are detected at 
similar luminosities; the remaining 6 have not yet been ob- 
served by Chandra or XMM-Newton (apart from snapshots, 
too short to provide meaningful constraints). 

For our fitted luminosity of SN 1941C, we infer a mass- 
loss rate from the progenitor star « 5 x 1O~^M0 yr~^ 
at « 55,000 yr before the SN (Soria & Perna 2008). For 
SN 1968D, the inferred mass-loss rate is ^ 8 x 10~^Mq yr"^ 
at « 30, 000 yr before the SN. The X-ray luminosity and 
mass-loss rate of SN 1968D suggest that this SN is evolv- 
ing to a remnant similar to Gas A. The similarity between 
SN 1968D and Gas A was also noted from radio observations 
(Hyman et al. 1995). The radio behaviour of SN 1968D was 
also found to be similar to that of another SN in NGG 6946, 
1980K (Hyman et al. 1995), which we have recovered in 
the Chandra observations from 2002-2004, at a luminosity 
4 X 10^'' erg s"\ in the 0.3-8 keV band (Soria & Perna 
2008). The presence of a power-law-like emission compo- 
nent above 2 keV in SN 1968D suggests that there may be 
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Figure 2. Fraction /^^^^ of Monte Carlo realizations of the SN 
sample for which the predicted pulsar luminosities are consistent 
with the observed limits of the corresponding SNe. This is shown 
for different distributions of the initial spin periods, described by 
Gaussians of mean Pq and dispersion crp^ . This sample includes 
only SNe older than 30 yrs. 

a contribution from a young pulsar and its wind nebula; its 
luminosity ~ lO^'^ erg s~^ is comparable to the luminosity 
of Grab- like systems (Possenti 2002), and in particular the 
Crab itself and PSRB0540-69 (Serafimovic et al. 2004). 



3 CONSTRAINTS ON BIRTH PARAMETERS 

We started from the results of a recent large-scale radio pul- 
sar survey (Arzoumanian et al. 2002): for a dipole approx- 
imation, the magnetic field has a log-normal distribution 
with a mean value log Bo(Gauss) — 12.35, asg ~ 0.4; the 
initial birth period distribution has a mean log Po(s) ~ —2.3, 
(jpg > 0.2. In order to test the resulting theoretical predic- 
tions for the pulsar distribution of X-ray luminosities, we 
performed Monte Garlo realizations of the compact object 
remnant population. The fraction of SNe that leave behind 
a NS has been theoretically estimated by Heger et al. (2003) 
as ~ 87% (subject to variability depending on stellar metal- 
licity and initial mass function). If an object is a BH, we 
assigned it a luminosity < 10^^ erg s~^. For the NSs, the 
birth period and magnetic field in our simulation were drawn 
from the Arzoumanian distribution, and were evolved to the 
current age. The corresponding X-ray luminosity was then 
drawn from a log-Gaussian distribution derived from the 
Possenti relation. When comparing the predicted and ob- 
served X-ray luminosity distributions, we need to consider 
the possibility of obscuration. However, we suggest (see dis- 
cussion in Perna et al. 2008) that neutral absorption does 
not strongly affect our results, particularly in the 2-10 keV 
band, because we considered only SNe older than 30 years 
(an analysis of all those older than 10 years gives the same 
statistical result); also, luminous X-ray pulsars would ionize 
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Table 1. Unambiguously- identified core-collapse SNc in nearby 
galaxies {d < 15 Mpc) until 1970, and luminosity or upper limits 
to the power-law component of their X-ray emission (upper limit 
to the luminosity of a possible young pulsaj"). All values derived 
from Chandra data, except where noted. 



SN ID 


Type 


Galaxy 


Distance 
(AllH-) 


ipo. 0.3-8 

(rrss >-) 


Notes 


1909 A 


II 


MlOl 


7.4 


< 5 X 10^® 




1917A 


II 


N6946 


5.5 


< 2 X 10^® 




1921B 


II 


N3184 


8.7 


< 10^^ 




1923A 


IIP 


M83 


4.5 


< 3 X 10^^ 




1937A 


IIP 


N4157 


15 


< 3 X 10^'' 


(a) 


1937F 


IIP 


N3184 


8.7 


< lO-'^ 




1940A 


IIL 


N5907 


13 


< 2 X 10^'' 


(a) 


1940B 


IIP 


N4725 


13 


< 3 X lO^'' 




1941A 


IIL 


N4559 


10 


< 2 X 10^^ 




1941C 


II 


N4136 


10 


5 X 10'''^ 


(6) 


1948B 


IIP 


N6946 


5.5 


< 2 X 10^® 




1951H 


II 


MlOl 


7.4 


< 5 X 10^® 




1954A 


lb 


N4214 


3.0 


< 5 X 10^^ 




1959D 


IIL 


N7331 


15 


~ a few lO^'' 


(6) 


1961U 


II 


N3938 


12 


? 




1961V 


Iln 


N1058 


9.1 


< 3 X lO^'' 




1962L 


Ic 


N1073 


15 


< 2 X 10^* 




1962M 


IIP 


N1313 


4.0 


< 2 X 10^6 




1964H 


II 


N7292 


15 


? 




1964L 


II 


N3938 


12 


? 




1966J 


lb 


N3198 


12 


? 




1968D 


II 


N6946 


5.5 


2 X 10^^ 


(b) 


1968L 


IIP 


M83 


4.5 


< 3 X 1036 




1969B 


IIP 


N3556 


14 


< 3 X lO^'' 




1969L 


IIP 


N1058 


9.1 


< 3 X lO^'' 




1970A 


II 


IC3476 


10 


? 




1970G 


IIL 


MlOl 


7.4 


^ 1037 


(c) 



Notes - ("'From XMM-Newton data; Soria & Perna (2008); 
Immlcr & Kuntz (2005). 



the entire mass of the ejecta on a time-scale between a few 
years and a few tens of years. 

We find (Fig. 1) that out of the 106 Monte Carlo real- 
izations of the sample, none of them predicts pulsar lumi- 
nosities compatible with the observed SN X-ray limits. The 
predicted high-luminosity pulsars, corresponding to those 
with the shortest initial periods, are not observed. This sug- 
gest that the mean initial period of the pulsar population 
is slower than the millisecond periods derived from some 
population-synthesis studies in the radio (Arzoumanian et 
al. 2002). A number of other recent investigations have 
reached similar conclusions. For example, the population- 
synthesis studies of Faucher-Gigucrc & Kaspi (2006) yielded 
a good fit to the data with a mean birth period of 0.3 s, 
a = 0.15 s, for logBo(Gauss) = 12.35. Similarly, the analysis 
by Ferrario & Wickramasinghe (2006) yielded a mean period 
of 0.23 s for a magnetic field of lO'^^ Gauss. We performed 
Monte Carlo simulations of the X-ray pulsar population us- 
ing those birth parameters and found them to be consistent 



with our SNe X-ray limits. We obtain a quantitative limit on 
the birth period by plotting (Fig. 2) the fraction of Monte 
Carlo simulations for which the luminosity of each pulsar 
is found below that of the corresponding SN (that is, the 
fraction of "acceptable" configurations), for different values 
of standard deviations. Mean periods shorter than « 40 ms 
arc ruled out at the 90% confidence limit, for any assumed 
value of the width of the period distribution (Perna et al. 
2008). This is also consistent with studies of young Galactic 
SNRs (Srinivasan et al. 1984). 

A possible reason for the lack of luminous X-ray pulsars 
is a turnover in the efficiency of conversion of rotational en- 
ergy into X-ray luminosity (relation between Edot and Lx), 
at high values of Edot- Another systematic effect that might 
affect our results is that a fraction of NSs may be born with 
a non-active magnetosphere. Similarly, the precise fraction 
of BHs versus NSs in the remnant population plays a role 
in our results. A larger fraction of BHs would alleviate our 
constraints on the initial spin periods, while a smaller frac- 
tion would make them tighter. Our results for the bulk of 
the pulsar population do not exclude that the subpopulation 
of magnetars could be born with very fast spins. But their 
spin-down times would be very short, and the spin-down X- 
ray luminosity would decline within the first few months, 
during which the SN is still too optically thick to let the 
pulsar X-ray photons go through. Therefore, our analysis 
cannot place meaningful constraints on this class of objects. 

Finally, our results have implications on the contribu- 
tion of young X-ray pulsars to the population of ultralu- 
minous X-ray sources observed in nearby galaxies. Earlier 
models predicted that a sizable fraction of that population 
could be made up of young, fast-rotating pulsars (Perna & 
Stella 2004). However, our study shows that the contribu- 
tion from this component, although expected from the tail 
of the population, cannot be significant. 
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